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The in-line force and the resulting drag and inertia 
coefficients for smooth and rough circular cylinders 
immersed in a sinusoidally oscillating flow at low Keulegan- 
Carpenter numbers (K) have been determined experimentally 
and compared with those obtained theoretically by Stokes and 
Wang. In addition, flow visualization experiments were 
carried out with oscillating cylinders in a water table and 
the stability of the flow was investigated. The results 
have shown that for very low values of K, the flow about the 
cylinder is laminar, attached, and stable and the drag 
coefficient is nearly identical to that predicted theore- 
tically. At a critical K, the flow becomes unstable to 
Taylor-Gortler vortices and the drag coefficient jumps to a 
higher value. Subsequently, the flow separates, becomes 
turbulent and results in a minimum drag coefficient. The 
subsequent increases in drag are attributed to vortex 
shedding. The inertia coefficient agrees with that obtained 
theoretically in the range where the flow in laminar. 
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I. 



INTRODUCTION 



The determination of the forces acting on bluff bodies 
immersed in time dependent flows has been of special 
interest in many branches of engineering. Special attention 
has been devoted to the wave forces on offshore structures. 
The extensive laboratory and field studies of the past 
decade have yielded extremely interesting results regarding 
the in-line and transverse forces acting on smooth and rough 
cylinders [Ref. 1]. However useful these investigations 
were, they were based on the use of an empirical equation 
(Morison's) [Ref. 2] and on the force coefficients obtained 
either through the Fourier Analysis or the method of least 
squares. In recent years, it became quite clear that 
further progress toward the understanding of the physics of 
the phenomenon must be based on analytical studies and the 
understanding of the kinematics of the flow. 

A theoretical analysis of the separated flow about a 
circular cylinder is an extremely difficult problem. 

Efforts in that direction through the use of the discrete 
vortex models did not lead to meaningful conclusions. It 
became clear that the key towards the understanding of flow 
separation, transition and vortex formation is in the very 
low Keulegan-Carpenter number (K) regime. It is because of 
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this reason that experiments were carried out in an 
oscillating water tunnel for force measurements and with an 
oscillating cylinder in a water table for flow visualiza- 
tion. The investigation resulted in findings beyond the 
original objectives. In the following, first the 
theoretical analysis and then the experimental apparatus and 

y 

measurements will be presented. These will be followed by a 
discussion of the results and conclusions. 
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ANALYSIS 



II . 



A. INTRODUCTION 

This section will present a condensed review of the 

theoretical studies regarding the oscillating flow about 

bluff bodies. First the flow about an infinite flat plate 

undergoing sinusoidal oscillations will be briefly 

discussed. Then Schl ichting ' s [Refs. 3 and 4], Stokes' 

[Ref. 5 and 6], and finally Wang's [Ref. 7] solution of a 

sinusoidally oscillating cylinder in a still fluid will be 

discussed. It should be emphasized that the case of a 

cylinder oscillating in a fluid otherwise at rest is 

identical to that of a fluid oscillating about a cylinder. 

The only difference comes from the 'buoyant' force acting on 

the cylinder in the direction of flow due to the pressure 

gradient needed to accelerate the flow, i.e., the inertia 

coefficient for the case of the oscillating flow is equal to 

(1 + C ) where C is the added mass coefficient for the 
a a 

oscillating cylinder, i.e., C = 1 + c: . 

IB a 

B. STOKES' OSCILLATING PLATE PROBLEM 

The flow about an infinite flat plate undergoing 
sinusoidal oscillations was treated by Stokes assuming that 
the flow is laminar and stable. With respect to the 
coordinate axis in Figure 1, the Nav ie r-Sto kes equations and 
the boundary conditions reduce to: 
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( 1 ) 



3 U _ 3 2 u 

3 t 3 y 2 

initial condition: u(y,0) = 0 

boundary conditions: u(0,t) = U Q cos ut 

u (® , t) = 0 

The solution of Eq. (1) via the separation of variables 
is given by: 



u(y,t) = U Q exp(-y/^ ) cos( u t - y ) (2) 

Letting n = y , one has 

u(y,t) = U Q exp(-n) cos(ut - n ) (3) 

"The velocity profile u(y,t) has the form of a damped 
harmonic oscillation, the amplitude of which is U Q ex P ) 

in which a fluid layer at a distance y has a phase lag with 
respect to the motion of the wall [Ref. 3]." Figure 2 shows 
the velocity distribution in the neighborhood of the 
oscillating wall. 

C. SCHLICHTING 'S OSCILLATING CYLINDER PROBLEM 

Schlichting [Refs. 3 and 4] considered the case of a 
cylinder oscillating at high frequency and small amplitude 
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in a fluid otherwise at rest and used the boundary layer 
approximation to the Navier-Stokes equations. This led to: 



a u 
3 t 



+ 





_ I IE + 

p 3 X 



2 



3 y 



which reduces to: 



I lE = iU . u 3U 

p 3 X 3 1 3 X 



(4) 



(5) 



for the outer flow, i.e., at the edge of the boundary layer. 
Combining Eqs . (4) and (5), one has: 





3 U 
3 t 



+ U 



3 U 

3 X 



( 6 ) 



This equation may be solved to the desired degree of 
approximation. The first order (linear) approximation 
yi elds : 



3 u ^ 



— V 



2 

3 U ^ 



3 U 

at 



(7) 



For the approximation to be valid, the convective 
acceleration must be considered smaller than the local 
acceleration, i.e., 
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3jJ 
3 t 



( 8 ) 
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3X 



Using the order of magnitude arguments, one has: 



2 

3 u U m 3u 

ax ~ D ' at 



m 



u 



It is seen that for the analysis to be valid, one must have 
U T 

m A 

— q— << 1 or ^ << 1 where U m is the maximum velocity in a 
cycle; T, the period of oscillation; D, the diameter of the 
cylinder; and A, the amplitude of oscillation. 

The solution of Eq. (7) with the following boundary 
conditions : 



at y = 0, u = 0 (no slip condition) 
at y = ®, u = U 

yields the following stream function: 

ip 1 (x,y,t) = U Q (x) f(n)e iu)t (9) 

from which one obtains the velocity as: 

u(x,y,t) = u o (x) [cos(ut) - exp(-n/ J~2) cos(ut - n/i/?) ] 

( 10 ) 
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which is seen to be identical to Stokes' solution for the 
oscillating infinite flat plate. 

A second approximation may be applied to Eq . (4) by 

retaining the convective terms. Then one has: 



3 u. 



3 t 



— v 



3 U 



3y 



„ 3U 8U 1 3U 1 

= u 17 - U 1 W - V 1 w 



(ID 



The stream function needed for the solution of the above 
equation is given by: 



* 2 (x f y f t) 




(x) 



< ^° 1 r c , x 2i ut . 

55“ u l f 2p <n)e + £ 



2s 



( n ) } 



( 12 ) 



where f 2 p represents the periodic motion and f 2s / the steady 
state. The use of the following boundary conditions on f 2 p 5 



at n = 0 u == 0 ==> f = f ' = 0 

at n = 00 f 1 = 0 



and on f2s : 

atn=0 u = v = 0 ==> f=f'=0 

at n = ® f ' = 0 
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leads to an important result. The boundary condition at 
n = ®r f' = 0, for the steady state component cannot be met. 
Instead a finite value is obtained: 



u 



2 




dUo 

dx 



(13) 



Periodically oscillating flow about a circular cylinder 
produces a steady streaming motion outside the boundary 
layer. In Figure 3 the streamlines about an oscillating 
cylinder in a still fluid are shown. This is for the case 
for a small Reynolds number with K -> 0 which leads to four 
recirculation cells and four stagnatiton points around the 
cylinder . 

The force acting on the cylinder has not been evaluated 
by Schlichting. 

D. STOKES' OSCILLATING CYLINDER PROBLEM 

Stokes, as one of his considerable achievements, 
investigated the flow about oscillating cylinders and 
spheres without resort to the boundary layer approximation 
[Refs. 5 and 6]. His solutions are then valid at an as yet 
unspecified small Reynolds number. His solution yields the 
following force per unit length of a cylinder: 






cos ut 




sin 



ut 

( 14 ) 
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which can be obtained by use of the asymptotic 



expansion of 



the Bessel function for large values of 
Stokes' solution will be taken up again 
discussion of a more recent solution. 



8 , ( 6 = D 2 / vT) . 

later following the 



E. WANG'S OSCILLATING CYLINDER PROBLEM 

Wang [Ref. 7] in 1968, investigated the flow about an 
oscillating circular cylinder by taking into consideration 
the effects of curvature which were neglected in 
Schl ichting ' s analysis. Using inner and outer expansion and 
the boundary layer approximation, Wang obtained the 
following force expression for the oscillating cylinder, 
valid to the order of S -2//2 : 



K 2 ° 



+ -1 - 
K 77 " B 



7T ^ 4 

•)COSUt - ~(1 + + - 

4 IT B / IT 6 / IT B IT B / IT B 



) sinut 



+ higher order terms 



(15) 



A comparison 
two terms in the 
Wang's solution 
third term. It 
for large B and 
since the third 
neg 1 ig ible . 



of Eqs . (14) and (15) shows that the first 

corresponding parentheses are identical and 
differs from that of Stokes only in the 
is noted that both solutions are valid only 
thus they yield virtually identical results, 
term in each parenthesis in Eq . (15) is 
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F. MORISON'S EQUATION AND THE STOKES AND WANG SOLUTIONS 
The normalized form of the Morison equation, with 
U = u m cos ut , [Ref. 2] may be written as: 



which was developed for the determination of wave forces on 
offshore structures and discussed in the literature in great 
detail during the past 30 years [Ref. 1]. 

The force expression obtained by Stokes and Wang may be 
expressed in the form of Morison' s equation by linearizing 
|cos e | cos e . 

From [Ref. 8], writing: 




2 

tt— C sin 9 
K m 



( 16 ) 
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(cos 9 I COS 9 



0 )cos 9 | cos 9 cos n9 d9 



n=0 




a + a, cos 9 + a 0 cos 29 + 
o 1 ^ 



with: a 



n 



0 n even 



n+1 



(- 1 ) 



2 



8 



n odd 



a 



n 



n ( n^ - 4 ) it 



_ 3 _ 
3 it 



19 



using a first term approximation: 



|C0S 9 | COS 9 = -r— COS 9 

J IT 



Then Eq. (15) may be reduced to: 



C = 2 + 4 7 = + A 7 = 
m / Ts b /Tb 



37T' 

C d - 2 BT 



i- + -A — 

it B 11 ^ 4nB/irB 
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case of 
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scillating flow 
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exper im 
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and 


rough 


cyl i 


nders . 
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Ill . 



EXPERIMENTAL EQUIPMENT AND PROCEDURES 



Force measurements were carried out in a large U-shaped 
water tunnel. It is 34 feet long, 25 feet high, and has a 
working cross section of 3 feet by 4.66 feet. Additional 
details of the design and operation of the water tunnel are 
given in [Refs. 9 and 10] . The flow in the tunnel is 
oscillated by means of a 1 HP fan and a butterfly valve 
system. The frequency of oscillation of the butterfly valve 
is adjusted by means of an electronic feedback control 
system, so as to oscillate the flow at its natural 
frequency. The amplitude of oscillations can be maintained 
at the desired level by opening or closing a gate at the 
exit of the fan. 

The test cylinders, approximately 3 feet in length, are 
mounted horizontally and supported at both ends by force 
transducers. These transducers can measure both the in-line 
and transverse forces. The electronic signals from the 
force transducers and from the amplitude gages are fed to 
two amplifiers and then digitized by means of an A/D 
converter. The digitized signal is then analyzed by means 
of a computer. 

The experiments for each cylinder were repeated at least 
twice to ascertain the repeatability of the data. 
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Two smooth cylinders with 6 = 1035 and 8 = 1380 and one 
rough cylinder with B= 1800 were tested. 

Flow visualization experiments were carried out in a 4 
foot by 8 foot water table. Circular cylinders of suitable 
diameter were oscillated by means of a slider-crank 
mechanism and a variable speed motor. The flow pattern 



about the cylinder was made visible by means of an electro- 
chemical technique [Ref. 11]. For this purpose a solder 
wire was inbedded along the cylinder length. Voltage of 
5-10 volts were applied between the solder and an 
electrode (placed away from the cylinder) to generate a 
metallic white smoke from the solder. A small amount of 
NaCl was added to the water to increase conductivity. 
Careful adjustment of the applied voltage and the repeated 
cleaning of the solder and cylinder surfaces were necessary 
for the successful visualization of the resulting flow 
patterns . 
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IV. PRESENTATION AND DISCUSSION OF RESULTS 

The drag and inertia coefficients were evaluated through 
the use of Fourier averaging [Ref. 1], as given by: 



Figures 4 and 5 show the drag and inertia coefficient, for a 

smooth cylinder, as a function of K for 8 = 1035 and 8 = 

1380. Similar data are shown in Figure 6 for a rough 

cylinder (k/D = 1/100) for 8 = 1800. Also shown in 

Figures 4 through 6 is the drag coefficient obtained 

theoretically by Wang, Eq. (15). 

Figures 4 and 5 show that C d follows the theoretical 

line for a K value smaller than a critical K and then 

c r 

increases almost abruptly to a higher value. Subsequently, 
C d decreases with increasing K and reaches a minimum at a K 
value of about 2. The variation of with K for K > 5 has 
been well documented in the literature [Ref. 3] and will not 
be dealt with here further. 





de 
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Figure 5 is similar to Figure 4 in many respects with 
the addition of a hysteresis in the drag coefficient. For 
very small values of K, C d followed the Wang-line and then 
rapidly jumped to a higher value. However, when the 
experiments were carried out by decreasing K at suitable 
intervals of K (with 15 minutes between each sampling) , C d 
remained at its higher value and then rapidly jumped back to 
the Wang-line at a K value smaller than that found in 
experiments with increasing K. The reason for the observed 
hysteresis is not easy to explain and it may be due to the 
precarious nature of the stability of the flow to the 
prevailing disturbances. It is tempting to think that in 
the experiments with decreasing K the level of the ambient 
turbulence is higher and this forces K cr to occur at a lower 
K. In the future, it may be desirable to wait longer 
periods of time between two successive changes of K. 

Figure 6 does not show any region where C d agrees with 
Wang's prediction. In fact, from the lowest K obtainable in 
the present experiments to the K at which C d goes through a 
minimum, Wang's prediction is considerably lower than that 
obtained experimentally. One will have to postulate that 
the flow has already become critical at K values smaller 
than about 0.45 and one is observing only the post-critical 
region of the drag coefficient. 

Figures 4 through 6 show that the inertia coefficient 
obtained experimentally agrees quite well with that given 
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by Wang for values of K smaller than K cr * It is noted that 
both the theoretical and experimental values of C are 



slightly larger than 2 due to the effects of viscosity in 
the boundary layer. It should also be noted that there is 
very little or no scatter in the C m data because the flow is 
in the inertia dominated regime. 

Figure 7 shows the RMS (root mean square) value of the 
normalized in-line force as a function of K for all the 
cylinders tested. The theoretical value of the said 
coefficient is given by: 



which is seen to agree well with the data for K smaller than 
about 8. This is due to the fact that the effect of on 
Cg(rms) is quite small and the value of (^(rms) is dictated 



The in-line force measurements and the resulting drag 
and inertia coefficients have revealed certain heretofore 
unknown facts which are both interesting and surprising. 
Interesting in the sense that the flow exhibits numerous 
transitions as K varies from about 0.4 to say 20. 
Surprising in the sense that one would not have expected 
transitions and hysteresis at relatively small K values 



m 



Cf ( rms) 




by C . 

m 
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particularly prior to the separation of flow which was 
expected to occur at a K value in the vicinity of 1.0 - 2.0. 
These two observations led to the flow visualization studies 
noted earlier. 

Figure 8 shows a sample photograph of the flow pattern 
about an oscillating cylinder. Clearly, Taylor-Gbrtler 
vortices form around the cylinder at regular intervals when 
the amplitude of oscillation exceeds a particular value. 
Figure 9 shows a plot of K versus e for smooth cylinders. 

For K values smaller than a critical K and for a given 6 the 
flow remains stable and attached. The symbol 'O' denotes 
the K value at which the Taylor-Gor tier instability has 
manifested itself. This instability gradually leads to 
separation and turbulence as K is increased. The symbol 'X' 
denotes the K values at which the flow has become turbulent 
in the boundary layer beyond which no clear vortices were 
observable . 

It is in light of Figure 9 that one can interpret the 
changes in C d in Figures 4 through 6. For 8 = 1035 the 
instability starts at K cr « 1.05, according to Figure 9. In 
Figure 4 the jump in C d occurs at a slightly lower K. 
Nevertheless, it is now clear that the rapid changes in C d 
in Figures 4 and 5 are a consequence of the instability 
observed in the flow. It is also seen that the K value at 
which the boundary layer becomes turbulent nearly 
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corresponds to the K value at which the drag coefficient 
goes through its minimum. 

The foregoing discussion has not addressed itself to the 
case of rough cylinders. Exploratory experiments with a 
sand roughened oscillating cylinder have shown that the 
effect of roughness is to reduce K nr > In fact, it was found 
that K cr - 0.4 at S = 1800. Such a low K value was not 
attainabile in the force measurements with the rough 
cylinder due to a number of reasons, the most important 
being the noise level in the signal. Thus, it is safe to 
assume that the flow for the rough cylinder has become 
unstable at a K value at about 0.4 or less and led to the 
rapid increase in prior to the lowest shown in Figure 6. 

In summary, it is now clear from the measurements and 
observations that the theoretical predictions of Stokes and 
Wang agree with the measurements as long as the flow is 
lmainar, attached, and stable. As the flow becomes unstable 
to Taylor-Gortler vortices, the drag coefficient is 
drastically affected primarily due to changes in the shear 
stress level and distribution over the cylinder. The next 
important change in the flow and in the force coefficients 
takes place when the flow separates and the boundary layer 
becomes turbulent. The inception of separation nearly 
corresponds to the occurrence of minimum drag. Beyond that 
the pressure drag becomes increasingly dominant and the 
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shear becomes negligible as K increases. Contrary to the 
previous assumptions, that Cj varies from 0 to its maximum 
value as K varies from 0 - 12, it is seen that there are 
both interesting and important phenomena in the low 
Keulegan-Carpenter number regime. The experiments described 
here, however difficult they may have been due to the 
problems associated with the measurements of very small 
forces, have shed considerable light on the low K flow 
reg ime . 
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V. CONCLUSIONS 



The investigation reported herein warranted the 
following conclusions: 

1. The theoretical solutions of the oscillating flow 
about a circular cylinder are valid only for attached, 
stable, and laminar flow; 

2. The flow becomes unstable to Taylor-Gbr tier vortices 
at a critical Keulegan-Car penter number for a given 
ratio of the Reynolds number to the Keulegan-Carpenter 
number, i.e., S; 

3. The drag coefficient shows a rapid jump at or near the 
critical Keulegan-Carpenter number. This is often 
accompanied by hysteresis; 

4. There is a second critical Keulegan-Carpenter number 
at which the boundary layer separates and becomes 
turbulent. This paricular value of K nearly 
corresponds to the K value at which the minimum drag 
occurs ; 

5. The theoretical value of the inertia coefficient 
agrees with that determined experimentally only in the 
attached, laminar, and stable flow regime; and 

6. Roughness precipitates transition and increases the 
drag coefficient. The effect of roughness on the 
inertia coefficient in the stable flow regime is 
neglig ible . 



29 



Figure 1. 



Coordinate Axis for an Oscillating Flat Plate 




Figure 2. 



Velocity Distribution for Stokes' Oscillating 
Flat Plate 
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Figure 3. Streamline Pattern of Steady Motion of an 
Oscillating Circular Cylinder 
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Figure 4. Drag and Inertia Coefficients versus K for 6 = 1035, Smooth Cylinder 
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Figure 5. Drag and Inertia Coefficients versus K for 6 = 1380, Smooth Cylinder 
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Figure 6. Drag and Inertia Coefficients versus K for B = 1800, Rough Cylinder 
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Figure 8. Photograph of Flow Visualization with a Smooth Cylinder 
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Figure 7. C f (rras) versus K for 0 = 1Q35, 0 = 1380, 0 = 1800 



* = transition to turbulence 
o= Inception of Instability 
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Figure 9. Critical Values of K versus g for Smooth Cylinder 
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